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Abstract. Large surface detectors might be sensitive 
not only with respect to extensive air showers induced 
by ultra high energy neutrinos but also to showers in- 
duced by hypothetical objects Uke microscopic black- 
holes. Microscopic black-holes might be produced 
in high energy particle collisions with the center 
of mass energies above the fundamental scale of 
gravity. These black-holes would decay rapidly by 
Hawking radiation into characteristic high multipUc- 
ity states of Standard Model particles and induce 
extensive air showers potentially detectable by a 
large surface neutrino detector. In this paper we 
study the possibiUty to detect microscopic black-holes 
exemplifying it in case of the surface detector of the 
Pierre Auger Observatory. The expected event rate 
is calculated for up-going and down-going showers 
induced by microscopic black-holes. Our calculations 
show a significant deviation of the expected rate 
compared to the rate expected by Standard Model 
predictions. The rates of up-going neutrinos are 
almost completely suppressed, whereas the rate of 
down-going neutrinos increase by a factor of about 
50 with respect to standard model calculations. The 
non observation of up-going neutrinos by the Pierre 
Auger Observatory in conjunction with a high rate 
of down-going neutrino-induced showers, would be 
a strong indication of physics beyond the Standard 
Model. 
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I. Introduction 

Searching for ultra high-energy (UHE) neutrinos 
(10^^ eV or above) emitted from astrophysical objects 
is one of the most challenging tasks in Astroparticle 
Physics. Neutrinos offer a unique opportunity to open 
a new observation window, since they are only weakly 
interacting and neutral. After having traveled cosmologi- 
cal distances without being perturbed and/or deflected in 
the interstellar medium, neutrinos behave as messengers 
of the most mysterious regions of astrophysical sources. 
Several theoretical models predict a significant flux 
of high-energy neutrinos as a result of the decay of 
charged pions, produced in interactions of UHE cosmic 
rays within the sources themselves or while propagating 
through background radiation fields. 

The expected neutrino rates possibly detected by large 
surface detectors do not only depend on the predicted 
flux of neutrinos but also on the neutrino-nucleon cross 
section. It has been noted f\\ that Earth-skimming z^^'s 



will generate upward going air showers when they inter- 
act in the earth crust. By contrast neutrinos of all flavours 
will generate deeply penetrating quasi-horizontal (down- 
going) air showers which are distinctive in having 
an electromagnetic component unlike hadron induced 
showers fT\. The rate for down-going v induced showers 
is proportional to neutrino-nucleon cross section, while 
the rate of Earth-skimming Vj- is not. Thus, the detection 
rate for down-going vs, and the rate for up-going v 
showers, react differently to variations of the neutrino 
cross section and tau energy loss. By comparing these 
rates one can therefore constrain significant deviations 
of neutrino interactions from SM predictions, see for 
example H for more details. 

The neutrino-nucleon cross section is related to parton 
densities in the yet unmeasured low Bjorken-x region 
X ^ 10^^. Some models even propose substantial 
modifications of neutrino interactions at high energies, 
including theories of TcV-gravity and production of 
microscopic black holes (BH), a domain that is to be 
tackled in the near future by the LHC. Measuring the 
flux of ultra high energy neutrinos would not only allow 
to put limits on cosmic ray production and propagation 
models. It would also probe fundamental interactions at 
energies that lie well above the TcV scale, and open a 
new window on possible physics beyond the Standard 
Model (SM) i). 

In this paper the scenario of microscopic BH produc- 
tion in neutrino nucleon collisions at the TeV scale is 
exemplified for the case of the surface detector (SD) 
of the Pierre Auger Observatory IS), which consists 
1600 water Cherenkov detectors with 1.5 km spacing. In 
addition to hadrons and photons the Auger Observatory 
is also sensitive to UHE neutrinos with energies above 
- IQi^eV 0, Q. 

The outline of the paper is as follows. In section II a 
short description of microscopic BH physics is given. In 
section III a full MC simulation chain starting with the 
injection of a neutrino into the Earth's atmosphere, its 
propagation, interactions and eventually the air shower 
production up to the actual response of the Auger SD 
array is described. Finally, in section IV results of the 
calculations event rates are presented. 

II. Black hole production 

In conventional 4-dimensional theories the Planck 
scale ~ 10^^ GeV (at which quantum effects of gravity 
become strong) is fundamental and the study of black 
holes lies beyond the realm of experimental particle 
physics. In theories where the existence of larger warped 
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extra dimensions is suggested, the 4-dimensional Planck 
scale is derived from the D-dimensional fundamental 
Planck scale, which can be of ©(TeV), bringing the 
possible observation of BH production and evaporation 
into reach of available experiments E), ID. 

Microscopic BHs might be produced in high energy 
particle collisions with the center of mass (CM) energies 
above the fundamental scale of gravity. When the impact 
parameter of two incident particles drops below the 
Schwarzschild radius of a BH with mass equal to 
their CM energy, BH formation should occur with a 
mass Mbh = V^- In case of a neutrino-nucleon collision 
the squared CM energy s is given by s = a; • 2m-[<sE^, 
where ton is the nucleon mass and Ei, the energy of the 
interacting neutrino. Within a semi-classical approach 
this suggests a geometric cross section of ct « Trr^ with 

1 

n + 2 

so that (7 cx s"+i, where n is the number of extra 
dimensions. 

The quantity AId denotes the lightest possible BH 
with a trans-Planckian Schwarzschild-radius. 

The total cross section for the formation of a BH in 
neutrino-nucleon collisions can then be written as 
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where Q) denotes the parton distributions functions 
(PDF) and is the minimal BH mass for which a 

semi-classical treatment of BH formation is expected to 
be valid and not understood effects of quantum gravity 
can be neglected. Due to the rapidly rising nucleon 
PDF at low X, the effective BH mass Mbh does not 
exceed several tens of TeV/c^, even for neutrinos with 
an energy up to 10^^ eV. 

Once produced, microscopic BHs are expected to 
decay within time scales of ~ 10^^^ s through three 
major phases: 

• The balding phase, in which the 'hair' (asymmetry 
and moments due to the violent production process) 
is shed. 

• The Hawking evaporation phase ifTOl . which con- 
sists of a short spin-down phase (the Kerr (rotating) 
BH loses its angular momentum) and then a longer 
Schwarzschild phase, which accounts for the great- 
est proportion of mass loss ifTTl . 

• A Planck phase at the end, when the BH mass or 
the Hawking temperature reach the Planck scale. 

Given the validity of the semiclassical description, a 
BH will mainly evaporate due to Hawking radiation and 
behave like a thermodynamical system with the temper- 
ature Th = . During the decay the Hawking temper- 
ature will rise as the BH mass drops. The lifetime can 



Approximately the 
decay may be treated as instantaneous on detector time 
scales, since for A/d 1 TeV and A/bh ^ 10 TcV, the 
BH lifetime t is smaller than 10~^^ s. During the decay 
process particles of all SM channels will be produced 
in a 'flavor-democratic' fashion with energies of order 
Th or above, typically at relatively high multiplicities 

III. Method 

In order to calculate the expected event rate from 
microscopic BH-induced showers at SD array of the 
Auger Observatory, a full MC simulation chain was set. 
Simulation consists of three phases: propagation and in- 
teraction of neutrinos inside the Earth and atmosphere to 
produce primaries able to initiate potentially detectable 
showers in the atmosphere; simulation of lateral profiles 
of shower developments in the atmosphere and, finally, 
simulation of detector response. 

The decay and the resulting particle spectrum of 
BH was simulated using a modified version of the 
CHARYBDIS code llT2t . The original version has been 
developed as an event generator for the production 
of microscopic BHs at the LHC. One characteristic 
feature to be pointed out is the inclusion of the recently 
calculated grey-body factors for BH production in extra 
dimensions llT3l . Grey body factors account for the 
fact that particles have to be transmitted through a 
curved space-time outside the horizon, and result in a 
modified emission spectrum from that one of a perfect 
thermal black body, even in 4 dimensions. For sake of 
simplicity and due to the fact that the decay balding 
and Planck phase are not well understood, only the 
Hawking evaporation phase for a non-rotating BH is 
implemented. In CHARYBDIS the fragmentation and 
hadronization process of the radiated particles is realized 
by an interface to the generator PYTHIA 1741 . 

There is no BH generator for collisions of neutrinos 
and nucleons available so far Among other generators 
also CHARYBDIS was originally designed for pp and 
pp collisions and motivated by collider experiments Uke 
the upcoming LHC. Modifications had to be introduced 
allowing for the study of neutrino-nucleon collisions. 
The neutrino has to be treated as a beam particle without 
substructure. This affects the calculation of the BH mass 
as well as the cross section, the handling of PDF and 
the initialization routine of PYHTIA. A Lorentz boost 
from the LAB system to the CM system before the 
collision and back after the decay was added to ensure 
the numerical stability even at highest energies. The 
implemented changes were verified by cross-checking 
the cross sections and distribution of the generated BH 
mass AIbh = \/xs (which depends on the involved 
PDF) calculated in our cross section generator, CSGEN, 
according to the cited literature against the results ob- 
tained with the modified version of CHARYBDIS. The 
cross-section generator, CSGEN, is a self-written tool to 
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Fig. 1. (Left panel) The BH formation cross section calculated for x^nin = 3. For comparison the SM CC u-,-N cross section 1151 is shown; 



(right panel) particle spectrum of a BH decay with Ei,^ 
the multiplicity (A/> = 614.4. 



10 eV, averaged over 100 events. The mean BH mass is (A/bh) = 6730.8 TeV, 



calculate cross sections. In CSGEN typical SM structure 
functions 111 5] were implemented as the parametrisation 
and other PDF included as the LHAPDF IHl FOR- 
TRAN libraries. 

All calculations were performed with the minimal 
black hole mass given by Xmin = ^^bh"/-^^d = 3 which 
corresponds to a lower cutoff where the semi-classical 
description of microscopic BHs is still valid, and the 
number of extra spacial dimensions set to 3 < n < 7. 

We determine the flux of neutrinos reaching the 
detector volume and initiating an extensive air shower 
(EAS). CSGEN is used to calculate cross sections, the 
distributions of involved kinematic variables and the tau 
energy loss for various interaction models and predic- 
tions ifTTl . The data are input to a modified version of 
ANIS ifTSl . |fT9l to simulate the propagation of incident 
neutrinos towards the detector and calculate the vertices 
of initiated air showers. 

The probability to detect neutrinos and microscopic 
BH by the SD array of the Auger Observatory is done 
by means of the packages PYTHIA and CHARYBDIS, 
which are used to simulate the production of secondary 
particles in neutrino-nucleon interactions that eventually 
initiate an EAS. These particles are input to the air 
shower simulation software AIRES ||20| to create shower 
profiles and footprints, which are then analyzed with 
the Auger Off line software framework 121)1 in order 
to determine the detector response and identification 
efficiency. 

IV. Results 

The actual BH formation cross section for neutrino 
nucleon collisions calculated with CSGEN is shown 
in Fig. [U (left panel) for different assumptions on the 
number of additional spacial dimensions ri. When the 
CM energy reaches values high enough to form a BH 
with the minimum mass Mg^" (which occurs between 
^ 10^5 eV and = 10^^ eV), the BH cross section 
rises rapidly, exceeding the SM cross sections by about 
two orders of magnitude at the highest energies. As 



an example the averaged spectrum of a BH decay at 
Ejj^ = 10^^ eV is shown in Fig. [T] (right panel). It is 
evident that the secondaries consist mainly of charged 
and neutral pions, and kaons. 

For the scenario of microscopic BH production, the 
neutrino interaction length in air is still larger than the 
atmospheric depth, so that BH showers can be initiated 
deep in the atmosphere and hence can be distinguished 
from hadronic cosmic rays in the same way like neu- 
trino induced showers, i.e. looking for inclined young 
showers 0, ll22l . 

Neutrinos are able to penetrate deep into the atmo- 
sphere before interacting and generating a young shower 
close to the detector as opposed to old showers of 
hadronic or photon origin shortly after entering the 
atmosphere. At large zenith angles the purely electro- 
magnetic part of such old showers is usually absorbed 
within the first 2000 gcm^. Practically only the high 
energy muons reach the ground, especially in inclined 
showers. This results in a thin and flat shower front 
which generates a short detector signal, lasting only 
a few ten nanoseconds. Young showers however reach 
the ground with a significant electromagnetic component 
still existent, showing a curved and thick shower front 
at ground that leads to broad signals with durations 
of up to a few microseconds. Together with the time 
information of the particles detected at ground and the 
elongated shape of the footprint, young inclined showers 
can be identified and their origin eventually attributed 
to neutrinos. The larger the considered zenith angle, the 
more pronounced are these features. 

To calculate the expected event rate from microscopic 
BHs we have defined a set of cuts where the number 
of events passing is maximal while the background 
contamination (due to hadron-induced showers) is kept 
minimal. A method similar to the one presented in ||6l 
was used. Applying the cuts to simulated neutrino and 
microscopic BH showers yield to the neutrino/BH iden- 
tification efficiency, which is defined as the ratio of 
showers triggering the detector and passing the cuts over 
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Fig. 2. The SD identification efficiency maps calculated for down- 
going i/c neutrinos producing microscopic BH. The efficiency is ploted 
as the function of shower energy Es and the injection height h. 

the total number of simulated AIRES showers. In Fig. |2] 
an example of identification efficiencies is shown. 

Once the identification is known the event rate can 
be calculated. The total observable rates (number of 
expected events) were calculated according to = 
AT X J^^""" A{E^) X <!>{E^) x dE, where is the 

isotropic j^-flux, AT the observation time and A{E^) 
the acceptance for a given initial neutrino energy, E^. 

In Tab. Uthe rates (number of events per year), for dif- 
ferent injected z^-fluxes are listed. The rates labeled with 
"WB" are obtained for the Waxman-Bahcall bound 1231 . 
Other rates are calculated for the GZK flux ll24l and a 
flux due to Topological Defects (TD) ||25l . In addition in 
Tab. U the ratio R between expected up- and down-going 
neutrino events is presented for different interaction 
models. One finds that this ratio is indeed sensitive to 
the underlying neutrino-nucleon cross section. A higher 
cross section results in a decrease of up-going and 
increase of down-going event rates, and vice versa. But 
given the fact that for a non-exotic interaction model 
more than ten years of data-taking (depending on the 
incoming neutrino flux) might be necessary to detect 
a single down-going neutrino, the statistical relevance 
of the measured data will at best allow to put a limit 

TABLE I 

The expected event rates and the ratio between up- and 
down-going event rates for different neutrino flux and 
cross sections. 



Scenario 


flux $ 


rate 1/yr 


rate 1/yr 








up-g. 


down-g. 


ratio R 


GRV92nlo 


WB 


0.44 


0.22 


2.0 




GZK 


0.12 


0.05 


2.5 




TD 


0.49 


0.32 


1.6 


HP 


WB 


0.23 


0.26 


0.9 




GZK 


0.07 


0.05 


1.2 




TD 


0.25 


0.39 


0.7 


ASW 


WB 


0.54 


0.13 


4.2 




GZK 


0.15 


0.03 


4.8 




TD 


0.61 


0.17 


3.6 


BH, n=5 


WB 


0.02 


11.30 


^ i 

~ 5^0 
~ 3^ 
~ 1 (Kin 




GZK 


0.01 


2.01 




TD 


0.02 


19.31 



on the cross section or energy loss model. In the case 
of BH production the picture looks different: the non- 
observation of up-going neutrinos in conjunction with a 
high rate of down-going inclined air showers, initiated 
deep in the atmosphere, would be a strong indicator of 
physics beyond the SM. 

V. Conclusions 

A complete MC simulation chain to study the mi- 
croscopic BH-induced showers has been presented. Our 
calculations show a significant deviation of the expected 
event rate in comparison to the rate calculated by SM 
predictions. The non observation of up-going neutrinos 
by the Pierre Auger Observatory in conjunction with 
a high rate of down-going neutrino-induced showers, 
would be a strong indication of physics beyond the SM. 
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